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THE DESIGN, ANALYSIS , AND EXPERIMENTAL 
EVALUATION OF AN ELASTIC MODEL WTNO 


Dy 

Ralph K. Gavin, III, and Chavalit Thisayakorn 
(I ) Introduction 

It is common practice in preliminary static aeroelastic analyses to 
estimate elastic increments in stability derivatives by utilizing a 

clamped- vehicle stiffness matrix in conjunction with an aerodynamic influence 

✓ 

matrix determined for the configuration by use of the linear, inviscid 
aerodynamic theory. These aeroelastic estimates are normally refined by 
’freeing' the structure and including vehicle mass effects based upon the use 
of mean axis vehicle coordinate systems, [l] The computer program FLEXSTAB, 
which was written hy The Boeing Airplane Company under the sponsorship of 
NASA-AMES, is a relatively sophisticated implementation of these basic ideas. 
FLEXSTAB has the capability of admitting two different types of structural 
representations. If the structural characteristics of the vehicle are 
adequately represented by the interconnection of beam elements, then FLEXSTAB 
accepts the beam El and GK characteristics and generates the required structural 
matrices. On the other hand, the structural matrices can be generated externally 
in a finite element program, without the restriction of "beam-like structural 
properties, and the results can be entered directly into FLEXSTAB. 

The primary objective of this study was to develop experimental data from 
a carefully controlled elastic model to be used in evaluating the effectiveness 
of aeroelasticity programs such as FLEXSTAB for vehicles of the orbiter class. 

In order to accomplish this objective at a minimum cost it was decided to 
utilize an existing rigid 5 % fuselage model for the 002 Orbiter configuration. 
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and to construct elastic wings for the model. The 002 Orbiter wings were 
straight with moderate aspect ratio and were therefore amenable to a beam-like 
structural representation. 

(II) Design and Fabrication of the Elastic Wing 

In view of the assumed beam-lilce structure of the 002 Orbiter wing, it was 
decided that the load carrying member of the model wing should be a beam with 
well defined. El and GK characteristics. It V'is further decided that the 
distribution of El and GK along the wing span should parallel the 002 ET and 
GK shapes provided to ASsM by the NASA-JSC Structures group. However, no 
attempt was made to scale the given (stiff) beam characteristics to the 5 % 
model. Rather the selection of El and GK was based upon the criterion that 
measurable deformations and stability derivative changes should occur under 
expected aerodynamic load. 

The basic design philosophy that evolved was that the aerodynamic loads 
should be transmitted to the beam via rib members. Further, in order to 
approximate the finite-element aerodynamic methods, the wing surface was 
segmented in a streamwise manner and each segment was rigidly attached to a 
corresponding rib element. Figure l is a planform view showing location of 
the Elastic Axis, as well as the basic aerodynamic sections for the wing. 

The MCA airfoil descriptions are given in Figure 2. As can be seen from 
Figure 1, the elastic axis is swept aft at an angle of 9*75 degrees from the 
vertical. 

Figure 3 depicts the basic construction of the elastic wing. Note that 
each panel section is made of low density Balsa Wood which is cemented to the 
supporting rib member. The region between adjacent panel sections is 
approximately 1/16" wide and is filled with an ultra soft Neoprene gasket that 
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iB cemented to the Balsa. Stress Calculations indicated that a high strength 
alloy of stainless steel could he used for the wing spar tut that heat treatment 
was required after machining. This process allowed a 2.5 safety factor on 
ultimate strength and about 2.3 on yield strength. (Ho allowance was made for 
dynamic loads.) Ribs were made of the same material as the spar for ease of 
attachment. The spar was designed to have a rectangular cross section so that 
the rib elements could be firmly attached to the spar. Table 1 contains the 
dimensional data for the spar. The computed values for EX and GK are also 
listed in this table. The cross sectional moment data which was computed by 
using the formulae [U ] 

O h 

K - bar - 6kg tanh(irb) . (l) 

3 „5 2a 

TT 

I = / z 2 dA = ba 3 (2) 

12 

The above symbols are defined in Figure h, 

(III) Design Calculations 

The purpose of this section is to describe the analysis methods that were 
utilized in the design of the elastic wing described in Section II. Two basic 
analytical tools were utilized in the evolution of structural specifications 
for the wing; namely the Doublet Lattice Aerodynamic lifting surface procedure 
and the finite element structural analysis method for beam members. In the 
following, we first discuss the analytical formulation for the general, 
time-dependent problem. After this broad notational framework has been 
established, the special static-aeroelastic and flutter problems are considered. 

A. The Structural Model 

As has already been pointed out, the principal load-carrying member in the 
wing structure is the spar element. The spar ib essentially a beam element with 



Spanwise Panel 
(location number) 

b - Thickness 
( inches ) 

a - Thickness 
( inches ) 

Lengths 

(inches) 

2-3 

0.391 

0.407 

1.078 

3-4 

0.369 

0.355 

2.828 

4-5 

0.273 

0.368 

2.828 


0.243 

0.310 

2.828 

: 6-7 

0.173 

0.300 

2.828 

7-6 

0.152 

0.246 

2.828 

8-9 

0.146 

0.188 

2.828 

9-10 

0.129 

0.152 

2.891 

Table 1. Spar Dimensions (See also Figure 3) 


Station Number of 
Spanwise Segments 

Values of El 
(lb - in 2 ) 

Values of GK 
(lb - in^) 

3 

63687.50 

40090.40 

4 

39837.50 

27166.72 

5 

32875.00 

15314.09 

6 

17575.00 

8730.85 

7 

11325.00 

3718.85 

8 

5468.75 

1990.35 

9 

2343.75 

1156.62 

10 

1093.75 

600.77 


Table 1. (continued) 

Values of El and GK for wing spar 

E = 2.9 x 10 T lb/in 2 

G = j ,12 x 10 ^ lb/in 2 


FIGURE 4. Definition of symbols for beam 

Cross - section. 



d| 1k . vertical 

displacement 

d 2l . -bending 
K rotation 

- torsional 
rotation 


FIGURE 5. Definition of structural nodal 

displacements. 
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rectangular sections of length 2.8 in. whose sectional area progressively decreases 
in the spanwise direction. There are 8 spar sections * each with a possibility 
of 6 physical degrees of freedom per section, implying a maximum of 1+'". 
structural degrees of freedom. However, in view of the planned testing of the 
wing at very small (less than ten degrees) angles of attack, it was decided that 
in-plane bending of the spar would he minimal and hence only torsional and 
normal bending degrees of freedom were retained for each beam element. The 
degrees of freedom associated with each element are defined in Figure 5. 

The elemental stiffness matrix for the member shown in Figure 5 therefore 
reduces to [2]. 


%K 


2 ,K 


l 3,K d i,K+l d 2,K+l d 3,K+l 


[K]. 


12 ei /r 
6 ei jt c 

X 


-12 EI / V 
; x 
i - 6 ei fi l 

x 


h El /i , 
x 


-6 ei x /£\ 
0 
0 


GK/Jt, 

0 , 
0 , 
GK/Jl, 


symmetric 


12 El /Z% 
x 

- 6 EI /l 2 , 

X 


k EI fl, 
x 


G£/Jl S 


( 3 ) 


A composite stiffness matrix can be generated by appropriately combining the element 
matrices in an underlying structural reference frame. In view of the fact that 
the spar is straight, the assembly task is quite straight-forward in this case and 
can be accomplished by overlaying successive element matrices and adding 


overlapping terms, e.g. , 
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6 

A 


6 

Of course the left end of the inhoard spar element is constrained to have no 
degrees of freedom since the flexible wing is attached to the fuselage at this 
point , 

A generalized mass matrix can be formulated for the flexible wing by using 
the results given in [2], i.e. , the mass matrix for the element in Figure 5 is 






d 3,K d l,K+l d 2,K+l d 3,K+l 


[M] 


K 


13 

35 


lit 

210 


5A1 2 

I 

x 

10A1 


r 

105 


21 

x 

15A 


3A 


symmetric 


_9 

* 

H 

VO 

131 

i 

X 

0 

13 + 

61 

X 


t 

70 " 

5A1 2 

W 

10AC 

35 + 

2 

5AC 


» 

-131 . 

X x 

-l 2 

I 

X 

n 

-ill 

I 

X 

,2 21 
1 X 


m + 

10A1 

TW 

30A 

U 

210 “ 

10A1 

105 15A 


0 


0 


6a 

0 


0 

3A J 


(b) 


As in the case of the stiffhess matrix, assembly of a composite mass matrix for 
the wing can be accomplished by appropriately combining the matrices, 

K = 1,..,8. The effects of rib mass were included by simply adding lumped mass 
and inertia terms to the diagonal elements of the spar mass matrix. A summary of 
rib mass properties is given in Table 2. 
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Station number 
of ribs 

Distances of ribs’ 
e.g. off from the 
elastic axis (in.) 

Hib mass _ 

(slug) x 10 

Polar Mass 
moment of inertia 
(slug- in 2 ) x 10~ 2 

3 

0.706527 

0.243999 

0. 140291 

4 

0. 651691 

0.201206 

0.098401 

5 

O.596856 

0.163412 

0.067019 

6 

0.542019 

0.130345 

0. 044076 

7 

0.487183 

0.101735 

0.027786 

8 

0.432345 

0.077311 

0.016626 

9 

0.377509 

0.0568 

0.009311 

10 

0.322672 

0.039931 

0.004781 


Table 2. 


Data on steel ribs of the wing 



Since the wing spar is swept at an angl'-, A, of 9*75 degrees, it is necessary 
to transform the spar structural matrices into a coordinate frame compatible with 


1.1 m * .mic 1* 

rame 

* 

Thir: Is <■ 

n.;:T ly accomp I 

11 short 

by 

performing the following 

trans rormati on 








[K], 

= [T] 

[K] [T]' 1 ' 






! 

r 1 

0 

0 0 

0 


0 ■ 1 

| 




1 

i° 

cos 0 

-sin 0 0 

0 


0 i 




0 

sin 0 

cos 0 0 

0 


0 


where T 

I 

l 

1 ° 

0 

0 1 

0 


0 



| 

1 0 
I 

0 

0 0 

cos 

0 

-sin 0 



I 

!_o 

0 

0 0 

sin 

0 

COS 0 . 


B. Aerodynamics 

Development 






The basic procedure that has been used to treat aerodynamic forces is based 
on the Doublet -Lattice-Method (D.L.M. ) [3]. Fundamentally, the D.L.M. yields a 
set of Aerodynamic Influence Coefficients , [Dj, relating the assumed harmonic motion 
of the normal wash , {w}, at specified points on the wing surface to the pressure 
differential, {ACp},' across the wing. Specifically, the integral equation 

w(x,y ,z) = f f K(x-£» y-n, z-Si to *-(0 ACp da (5) 

lifting surface 

is approximately solved as 


{w> = [D] {ACp} 

(6) 

where the velocity normal to the oscillating surface is 

{W} * \ Be [{w> e lWt ], 
and the pressure differential is 

(7) 

{Ap} = Re [{ACp} e iuyt ]. 

(8) 

Let us now consider specifically the computation of the 

aerodynamic forces 

later use. 
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In general, ^ is the distance from the ,]th. structural node to the l/4 chord 
point for the ijth. panel and r. , is the distance from the jth, structural node 

x j ■— 

to the 3 A -hord point for the ijth. panel (See Figure 6). The total force 
acting on the ijth. panel is given by 


f . . - Ap. , A. . 
ij ij 1.3 


(9) 


and will be assumed to be acting at the center of the 3/4 chord point of the 
ijth . panel. In order to define the sense of the various forces and moments due 
to the aerodynamic forces, we must first define the assumed positive displacement 
at node j. This is done in Figure 7> the force tending to increase the d^ 
coordinate is 


f ij = l =1 * P ij A i J’ 


( 10 ) 


the force tending to increase d^ is zero and the force tending to increase d^ is 


p 

f 0 . = -I l. . Ap. . A.., (a moment) 
3j i=1 id id 


( 11 ) 


where we assume that l and r. . carry the sign of their x-coordinate location. 

i J ij 

We will further assume that there are Q spanwise panel rows and P ehordwise panel 
rows, implying that n» the dimension of {d} is 3Q. 

Let us now consider the computation of the normal wash W in terms of the 
displacements at the structural node points. By definition, the normal wash must 
be equal to the substantial derivative of the vertical displacement of the surface, 
(implying no fluid flow through the surface.) In particular, we are interested in 
satisfying this boundary condition at the 3/4 chord points for each panel. The 
vertical displacement at any point along the ehordwise centerline through node J 


is 


z d " d ld ‘ Xd 3d ’ 


( 12 ) 


where the definitions of Figure 7 have been used. The substantial derivative 


w j (t) = gV ‘ a « - * d 3J - H ~ d 3r 


(13) 


is therefore 
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FIGURE 6. Span-wise slice of wing showing 

definition of symbols. 



J 3 


FIGURE 7. Definition of positive sense for 
displacements at node j. 




Consequently, the downwauh at the 3/4 chord or the i.jth. panel can "be written. 
(Assuming no chordwise deformation). 


W (t) = d, - r l - U.d-., 
i.1 ij ij 3 3 3d 


Further, if is assumed to bo oT exponent] n. L form. 


W. ,(fc) = W e 
1 ,1 A ,1 


iwt 


(l4) becomes 


W 


ij = d iJ (1 " ,) ' r ij d 3j (l “ ) ‘ D « 4 


3j 


(14) 


(15) 


Let us now return to (6) and establish the mechanism for computing the 
generalized nodal forces. Denote 


{w> = {w 1:L W^ 2 .... I V 21 W 22 .... w, 


2Q 


W P1 W P2 ' * * * W PQ } 


T 


(16) 


In (6), [D] is therefore a matrix of dimension PQ x PX. Let the force vector 
{f}, the pressure vector, {Ap}, etc. be defined using the same ordering as (l6). 
In addition, let the generalized forces at node 3s. be ordered as 


f,, = vertical deformation force 
Ik 

f = bending deformation moment 

&iv 

f*. = torsional deformation moment, 
3k 


and define 


{f} - {f y] f 21 f 31 


f f f 

12 22 32 


p p p . \ 

r lQ 2Q 3Q J * 


i.e. , if} denotes the vector of forces acting at the structural nodes. We must 
S' 

now develop appropriate transformation matrices so that {f can be calculated 
in terms of nodal displacements (d). From (10) and (ll), we can write 
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where the matrix [g] is defined as 



where [H] and [E] are defined hy 



i 



XU 



If we now combine (18), (20) and (2l), we can finally write 

{f> = -2= [G] [D(w)] _1 {[H] + ico[E]} {d} 

00 


( 2k) 


C. The Static Aeroelasticity Problem 

The bulk of the analytical work conducted during this study involved the 
estimation of elastic deformations under steady flow conditions. In this case, 
to is set to zero in the aerodynamic matrices and D is defined over the field 
of real numbers. Hence (2.4) reduces to 

{f} = [G] [D]" 1 [H] {d} (25) 

S ^ 

00 

Gravitational loading of the wing was neglected because it was found that these 
forces were small compared to expected aerodynamic forces. The result of 
combining (25) with the composite stiffness matrix whose development was 
outlined in Section III-A is the following expression 



XI 


[K] {d} = ~ [G] [Dr 1 [H] ((d) + a {1} ) {26) 

CO 

where (1} is a vector of length 3Q each of whose components is a one. In 
effect, the term in parentheses on the right hand side of { 26 ) contains a term 
dependent upon the elastic deformation {&} and a uerm dependent upon the rigid 
angle of attack, a. Now it is a simple manner to indicate the solution to (26), 

e.g. 

{3} = ([K] -ip [0] [D]- 1 [h]) _1 ip fa] [D]" 1 [H] {1} a (27) 

CO DO 

A computer program was written to implement (27) using the planar doublet 
lattice procedure (vortex lattice for steady (Tow) to calculate the deformed 
wing shape under a specified Mach number, q y -ind rigid angle of attack, a. All 
theoretical aeroelasticity results described in Section IV of this report, were 
obtained via this procedure. A listing for this computer program is given in 
Appendix A. 

D. Flutter Problem 

The basic problem in flutter analysis is that of determining if the wing 
will develop oscillatory motions under test conditions , This requires the 
inclusion of appropriate mass matrices and an unsteady aerodynamics capability 
into the existing computer code). 

By utilization of matrix structural analysis methods, discrete mass and 
stiffness matrices have been developed for the wing under consideration 
(See Section III). The resulting differential equation assumes the classical 
form 

oo 

[M] {d(t)> + [K] (d(t)> = {F A (t)}, (28) 

where [M] = mass matrix 

[K] = stiffness matrix 



{ d { t ) > = displacement vector for the structure 
{F^(t ) } = the applied aerodynamic forces . 

The basic approach was to ascertain those frequencies at which harmonic motion 
can exist as a solution to (28). 


It is normally convenient in structural analysis applications to work with the 
frequencies and mode shapes associated wiLh the unexcited structural system. 

If we set to zero in (28) and assume 


{aft )} = fd> e i(Ilt , 

(29) 

then (28) becomes, upon rearranging: 

(m 2 I - [M] _1 [K]) {d} = 0 

(30) 


It is clear therefore that the natural frequencies for the free system correspond 
to the eigenvalue of f Ml~"^ TK] and mode shapes {d> correspond to the eigenvectors 
of [M] fie] . Lst us denote the eigenvalues by au and the associated eigenvectors 

by {p^}, where i - 1, 2, ... , n. (n is the dimension of {d}.) Define 


[F] = [{p 1 > {p £ } {p^}] 

(31) 

2 2 2 

[ft] - diagonal oj^ .... ]. 

(32) 


The natural frequencies in radians per second are tabulated below: 


5.241 

X 

< 

1.305 

X 

10; 

3.121 

X 

io io 

7.982 

X 

lo! 

1.970 

X 

IDTn 

5.578 

X 

lo! 

1.320 

X 

10 *? 

4.788 

X 

10 

1.177 

X 

10“ 

2.383 

X 

10, 

7.810 

7.123 

X 

X 

10 9 

10® 

1.121 

5.220 

X 

X 

10! 

10, 

1*. 1*83 

X 

10' 

3.524 

X 

10, 

4.065 

X 

10® 

1.992 

X 

10 

1 

3.258 

X 

10 9 

5.589 

X 

!0 

2.151 

X 

10® 

1.059 

X 

10 

1.789 

X 

10® 

9.588 

X 

10' 


* 

Table - List of squared natural frequencies - ( rad/sec )‘ 
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If we return to (28) and assume that the solution (d(t)} is written as a linear 
combination of the basic vectors of [r] with time varying coefficient, we have 

n 

{a(t)} * E <j> (t) {pj (33) 

i=l 

The $^(t) are unknown scalar functions of time. We can place (33) into vector- 
matrix form by defining 

{<j>^(t)} = {^(t) • tj> 2 (t) <i> n (t)} (3*0 

and writing 

{d(t)> - [P] { 4> (t ) > (35) 

The substitution of (35) into (28), with some manipulations, yields 

+ [Q] {cj)(t ) } = [P]" 1 CM]” 1 {P A (t)}. (36) 


It follows therefore that by the introduction of the coordinates ijo , i = 1, 2, 
n, the left hand side of (36) is decoupled and hence amenable to straight 
forward solution. However, the right hand side of (36), which has only been 
written in functional form to this point , is in fact a rather* involved function 
of {d(t)}. With the substitution of (2b) into (36), the dynamical equation for 
the elastic wing finally becomes 


where 


f-m 2 I + [Q]} HU)} * {[S(i.i)] + iw[T(u>)]} f<|>(w)} (37) 

[S(u») ] [P]- 1 [Ml' 1 [G] [DU)]' 1 [H] [P] (38) 

CCj 

[T(m)] = [P]” 1 [m ] -1 [G] [D(o ))]" 1 [E] [P]. (39) 

OD 


The inverse Fourier Transform of* (37) is 


U(t)} - [T(t)] * U(t)i + [fi] C«f*(t ) } - ES(t)] * U(t)} = 0 (bo) 

where 

ri3(t>] - 1 1 !i(oi) n 

ITU)] * ; ' J 1 :t (m> J ’ 
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and * denotes the convolution operation. The investigation for instability or 
sustained oscillation inodes for the elastic wing is now reduced to that of 
finding the location of the roots for the characteristic equation of system (^0). 
Unfortunately, closed form expressions are not available for the elements of 
[ T ( t ) 1 and [S(t)3 due to the fact that [D(m)] can only be computed numerically 
for specified values of oi. Tests were made to determine if any of the lower 
natural frequencies corresponded to eigenvalues for (40) and it was determined 
that they did not. 

( IV ) Description of Tests 

Essentially two types of tests were conducted on the elastic wing after it 
had been fabricated at the Texas A&M Model Shop. 

(A) Static Tests 

The vortex lattice program was utilized to compute the aerodynamic loads 
that could be expected to act on each streamwise row of panels if the wing 
was rigid. The loads were computed for an angle of attack of 5° and 

« 80. lbs/ft. Weights equal to these spanwise loads were attached to the 
l/h chord point for each panel section. Then vertical displacement measurements 
were made for each rib at 25% and 75% of rib length. By using these measurements 
and data on the unloaded wing position it was possible to calculate the wing 
elastic pitching rotation. 

(B) Wind Tunnel Tests 

A series of six wind tunnel tests were conducted at the Texas A&M University 

Wind Tunnel for two orbit er configurations. The conditions were: 

q = 50 lb/ftp ; a = 2°, 5°, & 8° 
q m - 80 Ib/ft ; a = 2°, 5°, & 8° 

The elevon setting was held at zero for all tests. These tests were conducted 
for the Orbiter with elastic wings and they were repeated for the Orbiter with 
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identical rigid wings. The standard force and moment data, C^, C^, and C^, were 
recorded for each sequence of tests. In addition, a Cathatometer was utilised to 
make vertical displacement measurements at each spanwise rib location for both the 
leading and trailing edges of the wing. The Cathatometer was instrumented with a 
potentiometer so that displacement readings could be automatically read into the 
wind tunnel digital computer. One difficulty that was experienced during the 
conduct of the test sequence was that the Cathatometer had to be moved in order to 
make both leading edge and trailing edge measurements. Due to the unevenness of 
the floor, the reference point was therefore shifted causing some difficulty in 
checking test repeatability. 

(V) Comparison of Analytical and Experimental Results 

As indicated in Section IV, a series of static load tests were conducted to 
verify that the structural model used for the wing in the computer program was in 
good agreement with the elastic deformations actually given by the wing. Figures 
8 and 9 show curves of vertical deflection and elastic twist about the y axis 
derived from both experimental and analytical procedures. It is clear that good 
correlation was obtained for both twist and displacement in the static case, 
implying that the mathematical characterisation of the wing was adequate. 

Figures 10, 11, and 12 summarize some of the deformation and force data 
collected during the wind tunnel program. Figure 10 reflects the expected result 
that increased rigid angle of at back yielded increased z - direction deformations. 
Further, for a given q^ and «, - deflections increased uniformly in going from 

the wing root to the wing tip. Figure 11 provides a comparison of for the 
elastic and rigid models obtained from experimental procedure. These curves are 
so close to rigid C obtained from the analytical procedure that the latter has 

iJ 

been omitted from the Figure. Finally, Figure 12 displays the elastic twist 

(increment in angle of attack) rib a = 8 degrees for q = 

2 2 

50 lb/ft and 80 lb/ft . It is interesting to observe that a small ir'vri.i v** 



TWIST ANGLE (TWIST) VERTICAL DISPLACEMENT (CM) 


3 4 5 6 7 8 9 

STATION ON WING 

FIGURE 8. Comparison of the static joad test 

linear deflections. Note that the 

deflection along the leading edge 
is in general larger than along the 

trailing edge. 



FIGURE9. Comparison of the static load twist 
angle. 


VERTICAL DEFLECTION (CM.) 
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STATION ON WING 

FIGURE 10. Plots of vertical deflections ( in 

centimeters) measured in the wind 
tunnel. 


ro 

U3 



Lift coefficient Cl 




FIGURE II ~to. Lift coefficient Ct versus angle of 

attack a, q =80 ( leading edge). 
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s 

angle of twist was measured at some of the inboard wing stations. Appendix B 
contains a tabulation of Wind Tunnel Test Results . 

Finally, we wish to offer a comparison of experimenpFj. and theoretical 
aeroelastic estimates. The theoretical method outlined in Section III-C was used 
to estimate static aeroelasticity effects. The reader is reminded that the 
fuselage section of the orbit er was represented aerodynamic ally as a flat plate 
whose streamwise length corresponded to that of the wing root. Figure 13 shows 
a comparison of estimated and measurements for the wing leading edge z deflection. 
This Figure shows that this particular set of data correlated quite well. 

However, from Figure 1^ we find that the theory estimates for elastic increment 
in angle of attack lie below those obtained experimentally. In an effort to 
determine if this difference was due to an inaccurate estimate of local center 
of pressure for each panel section, the centers of pressure were shifted > 
forward by 20% and a new solution determined. While this does give better 
agreement, it does not appear that, the solution to this discrepancy can be 
obtained by a simple center of pressure shift. 

(VI ) Discussion 

In this report we have described the design, fabrication, testing, and 
analysis of a quasi-elastic orbiter model. The elasticity properties were 
introduced by -constructing beam-like straight wings for the wind tunnel model. 

A standard influence coefficient mathematical model was used to estimate 
aeroelastic effects analytically . In general good agreement was obtained 
between the empirical and analytical estimates of the deformed shape. However 
in the static aeroelasticity case, we found that the physical wing exhibited 
less bending and more twist than was predicted by theory. Although the cause 
of this difference is yet unexplained, there are several factors that may have 


contributed to it: 



VERTICAL DISPLACEMENT (cm) 



3 4 5 6 7 8 9 10 


FIGURE 13. Comparison of the elastic deflections. 


ANGLE OF TORSION (degrees) 



FIG URE 14. Comparison of twist angles. 


29 


1) Inadequate aerodynamic representation. (See above discussion) 

2 ) Imprecise wind, tunnel measurements. Although the wing was relatively 
quiet during testing some induced vibration occurred. This vibration 
along with cathatometer operator error almost certainly induced an 
unknown measurement error. 

3) Structural Integrity. A continual problem that was experienced 
during testing was that the rib-spar weld joints could easily be 
destroyed by improper handling. This failure would explain spurious 
data points like the one at station 9 on Figure l4. 

4) In-plane bending. The elastic wing was designed to have 
approximately the same stiffness in-plane as normal to the plane. 
Since the linear aerodynamic theory provides an inadequate 
representation for drag, it was not possible to adequately 
model this effect. 

In summary, the resits obtained in this report indicate that the linear 
aerodynamic and theories provide an approximate estimate of aeroelastic effects. 
However, our results imply that the theory un derest imates (at least in this case) 
the elastic deformation in wing twist and hence the effects of elasticity on 
lift. We believe that further testing and more complete aerodynamic models are 
required to resolve this question. 
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APPENDIX A 

Aeroelastieity Computer Program Listing 


The following pages contain a listing of the computer program used to compute the 
elastic displacements of the model wing. 



L. 


//TB476 JOB (T907-.4-C — ,*10,001 ,C-> THISAYAKOPN » • JOB 958 

//*WATFIV 

C 

C ************** J^^KI**************************************** ■■*■ ********* ***** 

C ** DISPLACEMENT ON MODEL WING ** 

C ** CAUSED BY DYNAMIC LOAD ** 

C *********** *************************** ************************************ 

c 

C NOTATION 

C 

C PLSPAPt Jl*LfcNGTH QF J-TH SPAR (INCH! 

C MsMUMBCR CF S°ANWISF PANELS 

c NF=nuyr£B of .DEGREES of freedom fob each panel 

C EDP(JI a JTH. ELEMENT of THE ELASTIC DISPLACEMENT VECTOR 

C FDPIJJ =JTH. ELtMEM C c THE RIGID DISPLACEMENT VECTOR 

C 



IMPLICIT PFAL-P U-H,0*7> 

01 Mr NST CM E STI FF ( 6 , 6 » . f STI F H 6 ♦ 6 1 ,ST I F F { 24 ,24 ) , 0 IFF 1 24, 24) 
DIMENSION Vt'Cl (24> ,VEC2( 24) 

DI«I NS ION F.I IB ),GJ(B), EDP(24>,RDP(24) 

COMMON /A*- *17 RLSP4=m 
C OMMON f i RT-A 2/ r HE T A1 

COMMON /APE AS-/ QtNV( 60, 601 , Rl8UH( 10 ) ,R IBLHI 10) 
r C vw^ f i • . t 4 / PjCJ'jVI 48,601, AI24, 30) 


o 


40 

FORMAT 

(10*12.31 

to 


42 

FORMAT 

f 4 X, 8? 15.51 

11 


44 

fo'jmst 

C5FI5.6 > 

12 


46 

F ?K MAT 

lT?Ot I? » 

13 


48 

FORMAT 

CRT 15. 61 

14 


50 

FORMAT 

(8 -10. 6 ) 

15 


52 

f r»«/ T 

(® C T0.3) 

16 


54 

?n~ »( ~ 

(B c 15.3) 

IT 


70 

format 

< • 1 • » 

1 9 


7? 

F0kMA T 

(*-*» * 

19 

1 

74 

format 

( * — • , 4 X , 1 V A L U F S OF GJ I POUND-INCH SQUARED) 1 1 

20 


76 

FCF^AT 

(• -• ,4X.' Vit'JFS ~i F FI (POUND-INCH SQUARED ) 1 1 

21 


78 

FOP mat 

( '-*,4X,* LENGTH 0F EACH SPAR (INCH)') 

22 


80 

FCF.vr r 

(*-»,4X, 'ELASTIC DISPLACEMENT VEC T GR (ALONG Tv*? MAJOR A XI 



l 

> 

23 


B2 

fF.pmaT 

( »-• ,4K , 'HERE I S TMC STIFF-MATRIX') 

24 


B4 

FORMAT 

(•-* »T2U, 'REDUCED SIZE OF O-INVERSE MATRIX!) 

25 


96 

format 

(*-• ,T?0,‘ /-MATRIX' ) 

26 


99 

FORMAT 

('-' »T2<), 'MATRIX K- MATRIX A') 

21 


91 

FORMAT 

(7X,» INCH*. 12 3, 1 R AO I AN • , T 3 7* 'RADIAN M 

29 


94 

FORMA T 

<4X, 'ELASTIC DISPLACEMENT VECTOR (ALONG TH= LEADING EDGF) 

29 


95 

FORMA T 

(7X, 'CEHTT PETERS' ,T22, • DFGREE S • ,T37 , 'DEGREES • ) 

30 

C 

96 

FORMAT 

(4X, * 2LAST 1C DISPLACEMENT (ALONG THE TRAILING EDGE)*) 

31 


M=8 


32 



w P2=m+3 


33 



NF=fe 


34 



M? = {<'!F/2» 4 ' 

35 

r 


THETtl= 

0.16972 

* 

36 

Lr 

(* 


CALL OL 

m (°E7A, SO, ALPHA, N) 


c 


REDUCE 

INVERSE D- MATRIX TH THE WORKABLE SIZE 

37 





36 



MNT=( m + 2)*sj 



U) 

ro 


' 39 



U=0 

40 



12=0 

41 



DO 890 J = 1 »N 

42 



12=12+2 

A3 



DO 890 1=1 . r 

44 



I J-l.l + 1 

AS 



12=12+1 

46 



KL-0 

A 7 



00 888 L = l,N 

AS 



DO 888 K=1,MP2 

49 



KL=K l +1 

SO 



SEOlNVt 1J.KL »=0!NVtI2,KL J 

SI 


888 

CONTINUE 

52 


90Q 

CONTINUE 

S3 



PRINT 94 

SA 



00 897 1=1, KN 

55 



PRINT 46, I 

SA 



PRINT 42, IP EPiNVII, JI,J=1,MNT> 

57 


897 

CONTINUE 

59 

L 

c 


CALL A« AT (M,N,Pi3,;LPH4tS0, 8ETA,HSTAPJ 


c 


PRINTOUT T Ht : A-MATRIX 

59 



PRINT 88 

60 



on ^oo i=i, «3 

61 



PRINT A6, I 

62 



PRINT 40, (All , J) , J=l, MSTA° 1 

63 

c 

900 

CONTINUE 


L 

C 


READ IN IMPUT-DATA 

64 



READ 52, 

65 



PEAD 52, (FIU»*J = l,K» 

66 



PflO 50, (PLS D APIJ),J=l,W) 

67 



DO 7 1=1. M 

63 



gju i=ojmn.i2o • 

69 

i 

7 

CONTINUE 


l 

c 


PRINT OUT ! NPIJT-DATA 

70 



PRINT 70 

71 



PRINT 74 

72 



PRINT 54, (GJ( J) , J=!.M> 

73 



PRINT 76 

74 



PRINT 5 A, CEII Jl ,J*1,M1 

75 



PRINT 79 

76 



PF I NT 48, tPLSFARCjl 

77 

U 


DO 1A 1=1, M3 

78 



DO 14 J=l,93 

79 


14 

ST IF C ( I ,J>=0. 

80 



mn1=m-i 


l# 

c 

- 

FORM GLOBAL STIFFNESS AND MASS MATRICES 

91 



CD 32 i:=i,"«i 

87 



I 3M3= I 1 +3-3 

33 



IFtll.GT.lIGD T n 519 

34 

L 


CALL STIFFS (SSTI c l.GJIl > , E I ( 1 », 1 , NF » 

35 



CALL TP FORM ("STTfil , THETA 1 ,NF ) 

46 



PALL STIFFS ( F ST T FF , GJC2 ),EIf?),2,MF! 

AT 



C.iU T-.FO?M. 1 F. ST TFF, THETA 1,NF) 



c 


ftft 


NF02-NF/2 

39 


DO ftftft I=1.NFP2 

90 


DP 6 ftft J=1,NF02 

91 


£STI*F f I.Jl=eSTlFF(I. J)*Ev . IF1U+NFD2»J+NF02) 

92 

646 

CONTINUE 

93 


GC TP 9 29 

9ft 

919 

rr'iiiNiJc 

95 

f* 

!1P1=U+1 

96 

l 

CALL STIFFS ISSTIFF.GH UFD.FtUlPl t,HPl,NF) 

97 

r 

CALL TP FORM tESTIFF,THrTA 1 ,NF ) 

9 ft 

929 

r 

CCNT1NUS 

' ' 

l 

c 

COMBJNF ELEMENTAL MATRICES INTO GLOBAL MATRIX 

99 


OH 13 ll=I,NF 

100 


00 I- Kl=l,NF 

101 


STIFF 1 13M3+J1, I3M3»KU*ST1FF( I3H3+J1, I3M3+K 1 i +ESTIFH J 1, Kll 

102 

n 

CONTINUE 

103 

32 

c 

CONTINUE . ... - 


t 

C 

PRINT ocT TMF STIFF-MATRIX 

10 ft 


PRINT B.2 

105 


00 2 R 0 1 = 1 r'N3 

10 ft 


PRINT ft ft, I 

107 


PRINT 54, tSTIFFI IiJlt 

10 ft 

2 B 0 

r 

C 0 V T IN‘I C 


c. 

c 

c 

V A T R I X K- matrix a 

109 


0 " 310 I =1 ,^3 

110 


2J 310 J=L,M3 

111 


OIFF ( I t .H=STlf F(1»J1-AU,J1 % 

112 

310 

CONTINUE * 

113 

\ 

GOTO 320 

lift 


PRINT Bft 

•♦WARNING** EXECUTABLE STATEMENT FOLLOWS A TRANSFER 

115 


□ n 3U 1=0, M 3 

lift 


PRINT 6 6 , I 

1 17 


PRINT 5 ft, 4 01 F FC I , J 1 » J=l »M3 1 

lift 

311 

continuf 

119 

320 

C 

J* 

CONTINUE 

120 

CALL INVfT IOIFF, M3, 1.0061 


l* 

C 

G c NcP* L I ?FC 01 SPLACfMENTS 

121 


00 33? 1=1, M3 

122 


VfC?< !1 =0.0 

123 


i>9 332 J=l,MSTA9 

12 ft 


V-VCZt 1 1 =vsr.2t I l+AU, JlMALPHA 

125 

332 

CONTINUE 1 , 

12 ft 


f>T 3^-0 1=1. M3 

12 ? 


C C»I I 1=0.0 

123 


00 360 J=1,M3 

129 


cn»n ) = cnpf I 1«-DIFF| I » J 1* VEC 2C J) 

130 

360 

r- 

CONTINUE 


L 

c 

PRINT TUT 01 SPIACCNFNT VECTOR 

131 


°S t NT 7 2 


132 PRINT ft o 

133 p««1WT <>1 

134 PR.’-VT 44, lcD»tU ,L=l»«3> 

135 PRINT 72 

136 PRINT 96 

137 PRINT 96 

138 on 600 1 * 1 , H 

136 J=t!-t)+3+l 

160 K=II-ll*3+Z 

141 L=U-l»*3+3 

142 POINT 44, FD«( -11*2.540 +RIBUHU ) *ECiP ( L 1* ?. 540 ,Enp(KJ*57.2 

1 96780, EDPIL)*57.2 957 80 

♦EXTENSION* OTHER COMPILERS *AY NOT ALLOW EXPRESSIONS IN OUTPUT LISTS 

♦EXTENSION* OTHE* COMPILERS WAY NOT ALLOW EXPRESSIONS IN OUTPUT LISTS 

♦EXTENSION* OTHER COMPILER S V A Y NOT ALLOW EXPPESSTONS IN OUTPUT LISTS 

143 400 CONTINUE 

144 PRINT 72 

145 PRINT 96 

146 * PRINT 95 

147 DO 420 1=1, M 

14-5 J=CI-l)"3+l 

14« K*I 1-11*3+2 

150 L=( I— 11 *3+3 

151 PRIM* r, { , t pDPt ,JI*?--6 aO -RTRLHII I *EDPI L 1*2 . 540 ,FDPm*57.2 

1 95 50 ? C 0PU 1*57.295750 

♦FXTEMSinN^ OTHER COMP UPS *AY NOT ALLOW E v .-«c?SinNS IN OUTPUT LISTS 

♦EXTENSION* OTHER CO* 1 P ILF PS *AY NOT ALLOW EXPRESSIONS IN OUTPUT LISTS 

♦EXTENSION# OTHFr COMPILER'S WAY NOT ALLOW EXPRESSIONS IN OUTPUT LISTS 

152 420 CONTINUE 

C 

153 99 STOP 

154 END 

£ *,****•+»*. pii**,-***. Mir*' + ** + ** +**«***4 *****+:«•♦*•■♦**» ♦♦♦itr#-#****^**** 


155 SU u -Rnu7If'r D-L v (PCTfi , SO, ALPHA, N1 4 

C * 

C THE DOUBLET L A T T ICE PR CC FOURS 

C FOR STFADY-PLANAP, SUBSONIC, CriMPRf SS I PL ? PLOW. 

C 

C INCPTMEMAL n S Cl l Lf TrtR y POWHWASH FACTORS 

C C8Y PITTING TH r *U.UEL FUNCTION ^OR LIFTING FUNCTIONS WITH A PAR A3 OLA ) 

C 

C SYMBOLS : 

C -• L C HA =STATIC ANGLE OF ATTACH 

C r '*■' T L rl = L TC A T I Of- 1 Cr THE ?DOT POINT ON LS SUING SCGF 

C ° CCT"5=lCCAT II *1 CP THE ROOT POINT ON TRAILING EDGE 

C T I P LI' =L CCA TI O'-l G c THE TIP POINT ON LEADING cGGF 

C TIPIS =LCC/TIGU rr thf TIP POUT on TRAILING EDGE 

C M = N I J u ° p c f;F COLUMNS DC WING PANELS < ON ONE WING) 

C N =UU*ip--s ns PANELS PF= COLUMN 

C ( THF COLUMNS RIJN FROM THE LEADING TO TRAILING EOGF OF THF «NGI 

C fl =1 FMGTH C e ,S rM I-Wl MG SPAM 

C P=prc m =R*:PCFN'TAC.S uf chopowise location on the wing 

C PEi-CSf J J=PcOC£HTAGc OF SPAM-WISF LOCATION ON THE WING 

C MUl I , J) = LENGTH of CACH OOUOLCT LIN? 

C •.PLMI ,Jh*.®:J OF "aCH WING PANEL 

C XIII.JI =LCCAT ION H r EACH SENDING POINT OUT EACH PANEL 

C X ( T « J I ^LOCATION OF RFCSIVING PCTNT ON tACH PANEL 

C AX5LP = SLOPE OF MAJOR AXIS OF THF WING 

C XIAXdUXI COOCDINATFS OF EACH COU.UMN ALONG THE MAJOR- AXIS 



u> 

vn 


*0+MHO Nil 



156 

157 

158 

159 

160 
161 
162 

163 

164 

165 

165 

1.67 

168 

169 

170 

171 

172 

173 

174 

175 

176 

177 

178 

179 

180 
k.81 
18? 

183 

184 

185 


186 

187 

188 

189 

190 

191 
L92 
193 
196 

195 

196 

197 

198 

199 

200 
201 
202 


c ntr-LC^n »=p 9 Fssusf; coefficicnts cm fnf wing 

c va =5PFrn np sound 

C FMCH =cfFr SIP* AH FA C H NUMRFP 

c U INF =F«CF STREAM VELOCITY 

c 

IMPLICIT RLAL*P lA-M t u- 7 ) 

PTMFMSICN PCVSLP(6),StOPF(6),PFRCSl 11) ,PEPCC17* 

DIMENSION XI I» 16) ,XIPR{6),X IID(6 l ,X IPO I 6 1 • TANLAMI7 ) 

Dl vf =f!SI nN XI1C 71 ,xtn 17) , YdOl.ETACUO) ,x!fX(lGl .CCllIrtTfrllll 
fji^NsnN xi(6aa)tXic{6,io),x(6fio) 

P 1 Mr NSI CM OrLPl6,li)),OFLX(6,lO)t M6»101 

0t*ier-*S107i 0 ELC{ 6 »n) 

RE ft L*8 L f '6 n £ (6 I, LAMPAOI 6 )> MUt 6, 10), I LEFT, IR IGHT 
COMMON /AR r 73/ 0° SS < 60 , 401 » R IGUH 1 10) » R IRLH ( 10) 

C n,, .ttON ft*it 5/ , T-L.M6.» 10 )•*=,- (6,10) 

C 

200 Pf'n I 5*?10,FND=990) 800TLF *ROGTTF,TI PLE»T IPTE*8,.Nr ►.ALPHA 
210 FC-V6T (5 C1 , 0.3»’’I 10* FIO « 61 
PP I NT ? ,->0 

220 FORMAT 1*1*1 ‘ 

PRINT 220 

230 FORMAT I4*,*0AT4 INPUTS- *) 

POINT ? 40 f R r 0 T tF,?09TT6.TIPLFtTIPTF t P,N f N 
240 FQRwAT C4X, • RCr’lE 3 • , F 10. 4, 4X . • «OOTT £ = * , Flt).4,4X, *TlPtE= »,F10*4 
1 ,4X, ,T IP':= * ,= lQ.6 f 4X, * B = « ,F10.4,4X , *M = *,I4»4X.*N= *»I4) 

POINT 2 50 1 'L«*H* 

250 FORMAT ( , 4 X ,* ANGLE CF ATTACK IS *,F 10 . 6 > 

C 

VA=10°0.0 

W=0.0 

=0.0 

UTNF=25i).0 
FK'ifH=UINF/VA 
RFTt=0$CPT( 1.0-FMACH**? ) 

n*8»f»FTA . * 

* 90=0. 5t>f=0, 002503 867* ,| JI NF**2 

C 

PRINT 252, FM A C H » 8 FT A , SO 

252 FORMAT (T20»*«ACH NUMBER** ,F10.6,T60, * BFTA=* t F10 .6, T90, • DYNAMIC Pft 
lESSOPEf 0=* * FI 0.3 l 
C 

C PK"CALOULA T I OH Fr<0 The <~COPDlNAT.rS of THE sending and receiving elfments 
C CALCULATION OF THE *9 FA CF EACH WING PANEL 
IN=NU 
I M=M*1 

9 LAC 254 , l*>E»CS( J), J=l, IM| 

4F6P 254, tPF»CCm,I = l , I N) 

254 e'U'AT 18^10.61 
PRINT 2 SB 

259 FC^AT r *-• ,4X, ‘PERCENT CF SPANWISE PANFLING LOCATION* I. 

PRINT 305, (PEPCSIJl ,J=1,IM > 

PRINT 260 

260 C C' 5 .*'AT l *- , ,4X, 'PERCENT OF CHORDWISF PANELINO, LOCATION*-) 

PRINT 305, (OEPCCin ,I«1,1H) 

CO 270 .1=1,1 M 
FTf r ji = ?^p= s cs(j) 

270 CONTINUE 

DC 290 t=l,IN 

xii in = ip?cTTE-BOOTLEK-PFP.ccu i+aooTtr 


70.3 XI mi ) = ( TIP T*--T[PLE» *?•!<> CC ( Il+TIPLP 

2G4 TANLAM( I IMXnil l-XIIUl »/s 

205 2 SO CONTINUE 

206 C ( 11 =p-»nTT r -S10TI" 

207 C (If'I-TIf’Tr-TIPLf 

203 fX5LP = (TIf > ti-«r'1TLc+0«37 5O l MC(IMl-C( 1 ))!/« 

206 Axr'T : 5 = 5 7.?8 5 796*0ATAN( AXSLP ) 

210 PRINT 2 84 , AXOEG 

211 284 FQP«AT (*-«,*SLr lJ E HP MiJjP. AXIS =«,F16.6,* DFGREFS * ) 

C 

C SLOPE OP Liw t^tJPLFTS. SLOPE > WITH ITS ANGLE, LAM!) A • 

212 on 370 1=1 ,N 

213 DO 300 J-l t I* 

214 CUI=CMI-CCm-C'lMM*ppRCS( Ji 

215 CELCU, J» = (P'2CC( 1+1 l-P r PCClin*C(Jl 

216 300 CONTINUE 

217 XlIDU > =XIIU ) *0.75o *r>CLC»I,lt 

21 R XI00( I ) sXIC'l 1 1 ♦O.2'5 , O»0ELCt I « I M ) 

210 •%* SLC1PI-U » = (Xl^oi Jl-XllJI I IT/S 

220 L A'lnA II) =D AT4*J ( 5|. CPC 1 1 1 ) 

221 ItvoAOt !I =LAM0* (11*5 7.20 47 8 0 

277 XIlMIlsXH 111*0.750 *OFLCC I>1) 

223 ■ XlIlPH»s.<IE(! >*0.750 = i;?LCU# I'M! 

224 370 BCVSLP* ||=(XI0PI I »-XtIR(IH7« 

225 PPJNT 380 

226 3B0 FO°*AT ( , V* *• ANGL = GF EACH DOUBLET LlNf IIN DEGREES) * 1 

227 PRINT 305, I L AMD A D(T)«t s l«N) 

228 305 F0BMAT (4X,F16.6) 

C 

C 

c Nfj computing thp rfouireq coordinates 

220 IJ=0 

230 r>" 405 1 = 1, N 

231 DO 405 J=1,M 

232 MU( I*J1=(E T C(J + 1 I— FT A ( J ) 1 /DCOS ( LAMP A ( I I ) * 

233 \ INDEX I=2”(m-J) *1 

234 DELA1 I. i J=0 . 5DD1:=0ELC < I, m (FT A 1 .1+1 1-ETA( J I )*(C( 11*INDEXJ 

1 *C(IMI’ , (?'.J-1»»/(C(U+M1 

235 DFLAl I, J1=DFLAU .JI/3ETA 

23ft I J-=! J + l 

237 XKT,.])=XIIDU )*' Ts(.J )*SLOPF(I) 

238 CTACt JJ =0.5D*(rT.*.(j J+ETA(J+1> ) 

230 Yt J)*5TACIJ1 

240 X ( ! , Jl = XI I e ( I J *PC VSL P( 1 ) +YV J! 

241 PEL.X.f I, )l=O.50mr<rLCU ,J 1*CELC ( I , J+l) 1 

242 XICU . H-X( ! ,J J-).5U»0HLX(I ,J) 

243 XIAXI J) -O.E75U*-.' { U+AXSLP*Y(,J>*RDOTLE 

244 FL( I, JI=XtC(I. Ji-X!AX( J) 

245 Pit, Jt = XlI, JT-XUXUt 

246 a I 3 !!H { .11= I C ( I KC f J*l n / 2 . 0* 0.3 75 0 

247 F I'M.M j | -RI3UU(J) *0.62 50/0.3750 

248 408 C nt ! TIN J" 

2 40 PRINT 708 

250 308 FORMAT (•-*,' XI-LCCATIQNS ALONG THE MAJOR AXIS:-’) 

251 PPIUT 340, (XI AXt *> ,<=1.V1 

252 PRINT 310 

2 53 310 rOPV'T (•-» ,4X ,’HE= C IS THE VECTOR OF Mi- 

254 Of* 311 1=1, N 

255 PRINT 330 

256 *3 I NT 3 4J), (MU(I , JT,J=l, PI 



t|U ONWO* 


257 


25 * 

259 

260 
261 
262 

263 

264 

265 

266 

267 

268 

269 

270 

27 1 

272 

273 

274 

275 

276 

277 


276 

279 

280 
281 
282 
293 


284 


285 

286 
287 
28 3 

289 

290 

291 

292 

293 

254 


255 

296 


311 CONTINUE 

C POINT 314 

C 314 FORMAT { •-» , 4X, • X-LOCATI OMS OF EACH RECEIVING POINT*) 

C PRINT *"0 ’ 

C PRINT 340, |(X<I,j) , J*l. P) t 1=1 .Nl 

PRINT ?20 

320 FORMAT ( ' -• »4X» 'PANEL AP EA ON 0N« WINS') 

00 350 1=1, N 
PRINT 330 
330 FC>MAT I * — • 1 

PkJ MT 340, (O' 1 LA I T » J1 , 4= It H'l 
340 FORMtT 14X,30F12.6) 

350 CONTINUE 
PRINT 352 

352 FORMAT ( , T?o, 'Ll I , J ) WITH SIGN OF THEIR X-COORDI NATE LOCATION 9 
IE ENG C7PRIF0 ALONG* ) 

00 353 I=l,N 
pP I NT 3 30 

P e I NT 340, C Ft < I » J) , J=l , P) 

353 CONTINUE 
PRINT 354 

354 FORMAT T20, • R [I ,J ) WITH SIGN OF THFIP X-CDORDI NATE LOCATION 8 

IE TNG CARRICO ALONG' ) 

DO 355 1=1, N 

OF J NT 3 3& 

P D Tr|T 3 40, (®( I,J1,.J«1,M) 

355 CONTINUE 

C 

C NOW CALCULATE DR 5 { l > AND DRS<SJ MATRICES 
C PRINT 360 

C 360 FOR NAT ,18X,'DPSS* ,20X,'DRS1M 

Tl=1.0 

K L=0 . . 

DO 590 K=t ,N 

DO 690 L=1 ♦ M . * 

, KL=KL+l 

1 J=0 

C PRINT 415 

C 415 FORMAT I*-') 

C PRINT 363, KL 

DO 5R0 1=1, N 
C PRINT 315 

C 315 FORMAT 14X, » U 
DO 590 J=l,M 
IJ=TJ+l 

XMXI=X(K,L>-X!II.J) 

Y M ETA= Yf L ) - C TA t J ) 

YPF TA =Y 11. > +E TA I 1 1 
Xi'XtC=X IK,U-XIC< I.J) 

VMfTAC= YID-FTACUI 
YPETAC= Y(L ) + c T7C l J) 

ROOT- OSORTCMtH I, J )**?-?, O *( XMXI*OSIN<LAMOA< 1)1 
1 *YMETA*0C0SIL4MCM1 I) ) )*NU( ! , J » >X«XI* ^2+YMFTA **2) 

l RIGHT* MU ( I , J ) /( YMF'T.A* (YILl-cTAI J+l-l » ) M ROOT/ 1 V tL ) ~ ETA l J +1 1 ) 

1 - OSnPT t X'4X ! t M2*-Yf*FTfi+-2 J/YMFTA ) /I XMXI*OCQSI L AMDftl 1)1 

2 -YMfcTA’fCSTNTLAMDAtl) H ' U) 

RCOT= DSOP-i 1*1(1 , Jl-*2-?,0 *{XMXI*f)5tNCL6Mf)A<I )) 00 

1 -YPFfA*OCnS ( LA W DA( II > V*«Uf 1 , J ) +XMXI **2*YPE TA**2) 

ILEF T=CU1 I, J)/[YPrTA*( YtL)+ETAU+ll H- 1 ROOT/ I V ( L )«-£TA( J + M ) 

1 - OSORT f XM X 2 + YPETA **2 ) XYPE TA 1 /< XMXI40C0SC LAMDAl II). 


2*1 

258 


2 Q 9 

300 

301 
?02 

303 

304 

305 


306 

307 
303 
30 Y 

310 

311 

312 

313 


314 

315 
315 
317 

3 13 
315 
320 


32 L 
~22 
3?3 

324 

325 
3 *^ 
3 2 t 
’28 
3~5 
333 
3?1 
33 ? 
313 
3~4 
335 
3 36 

11 7 

315 


340 


34 *. 

34 ? 


2 ♦VPrTttnsiMLAMfUUl 1) 

C DPSS= ( I P I0MT + ! L' :p Tl’f'OCOS (LAMDA C I ) > 

C IF IW.FQ.O.OI GOTO 452 

OPSSIKL*IJ)*UPIGHT+ILeFTl«cr.OSCLANDACII I 
nr. S$(«L *IJ1=0BSSIKL,IJI*Q. 50*Df=LXri ,JI /4.0/3. 141552 653 60 
C 

C TFANSFnRM 3 AC* TO THE ORIGINAL Y-AXIS 

Y M £ 7A=YM£T£/ BJr T A 

YorT/.sYP£T*/8ETA 
Y M r T CC= Y^5rTAG/ p F ta 
YPETAC= YP cTAC/PET A 
580 CONTINUE 
550 CONTINUE 
C 

CALL JNVPT » nRSS , MN , l . 0 P 01 
C 

00 600 1=1, N 

on too j— 1 » m 

FL U,J)=FL M,J1/5FTA 
600 CONTIwiJF 
cnift ?oo 

4*0 CONTINUE 
RFTUPN 
=N0 

Q * = *5 * "t. *Jt i *•** ,M M If,,* B I? tt + ,e $♦ * + + **4 4 *<! 7 *K*Hl* * ^ 


S'JPPnurif|= S T ! FFE I = STI FF, E 5 J ,EPI ,N, NF ) 
ILLICIT RF*,L=R IA-H, 0 -ZI 
OINENSirf; 7 ST 1 F C (\F,NF) 

CQf|M~r’ UP. hi?./ PL SPAR (81 
C 

Pic =r L SPAR ( N 1 

FLF?=PLC**? 

RLC 2 = RLC **3 

C 

C epf?M ?c.vrf*,r; STIFFNESS MATRIX - 

' [ f: 40 I= 1 ,’jF 

40 J=t.NF 

rsTiPFf i, .n=o.o 

40 CrvTINilF 

FS’lFFt 1 , 1 1 = 1 2 , /R LC ? *£ E I 
r ?. T 5 F =1 1 , 21 =G./SLC 2 *FFI 
?S T I"( 1 , 4 l=-l 2 ./RLC 3 «-rfI 
SSTIfrl 1 . 5 »s 6 ./PLC?-r«r 
£ c TIFr(?,l)=FSTIF c <l ,21 
FSTIF?( ?, 2 l= 4 ./PLC^Fn 
TSTIFFt ?, 4 » = - 6 ./RLC?*lFI 
E S r T F F ( 2, 51 = 2 . /PLC*t’FI 
E S * T F F ( 4 . 1 )-CSTIFF(l t 41 
tS’IFFf 4 , 2 >-F:STIFP{?, 4 » 

= ST 7 ?F( 4 t 4 j=fTSTIFMl ,11 
F STIFF* 4 , 5 >=-A ./FLC?*FfI 
F 57 IFF( 5 , 1 }=CSTIFF( 1,51 
F STIFF! 5 - 21 =T STIFF ( 2 ,51 
ES T IFFf 5 , 4 )=r 5 TIF=( 4 , 5 I 
E«TTFF( 5 ,S) = FSTIPF (? f 2 1 

n 

c =rt3 M TnnSI'.TMAt. STIFFNESS MATPIX- 

= STIFF( 3» 31 =LG J/P LC 
"S T IFF(6,3 l=— EGj/RLC 




3«3 

344 

345 

346 


ESTIFFI 3.61 »E STIFF 16.3) 
f ST I FF I 6, 6 J = tSTI FFI 3 .3 » 

RETURN 

FND 

***************** *******:************** ******* ************* ******** 


347 

34R 

344 

350 

351 

352 

353 

354 

355 
3 56 
357 
36 3 - 
3E*i 
3PQ 
361 
36? 
36 3 
>6^> 

365 

366 
367 
36» 
564 

370 

371 
37? 

373 

374 

375 

376 


SUBROUTINE TRPGRM l A .THFTAT ,NF ) 

IMPLICIT REALMS lA-H.O-Z) 

DIMENSION A{NF t NF),TRANSC6 ,6 >,PR0D<6 .6 » 

C 

00 44 1 =1 »N F 
DP 44 J*1»NF 
44 TP ANSI I »J)=Q.O 
TPflNSn»n«1.0 
T«AN514»4)=1.0 
TP ANSI? ,?l=DCOS( THETA T) 

TPANS<2,3I=-DSIN( THETAT) 

IF CMS I 3 .2) =OSIM THETf.T ) 

TRANS 13,3 ) =OCP S( THFTAT J 
TPA*!S(5,6»sDC0SlTHETAt» 

TP*N5<5 i 61 =-OSIN( THE TAT) 

TFANS(6,5)= DSINITHETATI 
TF ANSI 5 ,6) = DtOSt THLT AT ) 

OP 51 I*1.NF 
OP 51 .1 = 1 *N C 
PFODII. J>=0. 
nn 52 k=1,NF 

52 pppDn,.n=ppeoti, j)+tRANscK t n*ACK,jj 

51 CGNTT MU" 

ro s? i=i .nf 
on 53 J =1 » NF 
A ( I * J) = 0» 0 
op 64 K = 1 »NF 

54 A II.J)=A ( I,J)+PPOCI I,K )*TRANS((t, Jl 

53 CONTINUE , • 

, r ftijp w 

Ef.O 

Q ************ ***->»**-:.*.* **y ********** ****** ************* ********** 


377 


SUBROUTINE AM AT I M, N.M3 , AL PHA , SO. FET/,MSTAR) 


C 


c 

SVMRCLS 

c 

THtTAl 

c 

th; taz 

c 

B 

c 

M-OTLE 

c 

PC'PTTE 

c 

tiplS 

c 

T I PTE 

c 

AATDFG 

c 

EXF I'Ll ) 

c 

F XM L2> 

c 

Ex F ( 1 3 1 

c 


c 

Y( .j) 

c 

X A X 1 J ) 

r 

XU)(JJ 

*K 

ARM (. J | 

c 

P.l C 

C ' 

Ff- Cl* Jl 


= ANGLE FROM ELASTIC AXIS ONTO THE Y-AXIS’ 

= ANGLO FPD« LINE OF LOADINGS ONTO THE Y-AXIS 
= SFMI-SPA*!MTS? LENGTH 

= L0CATin*1 OF T J-E RC5T POINT ON LEADING F^r.E 
= LPCA TJON OF THr FOOT POINT ON TRAILING EDGE 
= LOCUIv‘l OF THE TIP POINT (JN LEADING EDGE 

*inC4TID*: OF THE TIP POINT ON TRAILING EDGE 

=A f, GLE OF ATf.CK IN DEGREES 

^VERTICAL FOKCE LOAD ONTO THE J-TM. SECTION OF THE WING 

sHFNOING f<0 ,/c NT LOAD ONTO THt WING (EQUALS ZERO IN H£R£ ) 

= TOR SIGNAL MOMENT LOAD ONTO THE J-TH. SECTION OF THE WING 
ICAUSFO OY VERTICAL FORCE APPLIED OFF FROM THE MAJOR AXIS) 
-Y— LOCATION r-F EACH WING SECTION 
=X-LQCATICN ALONG THE MAJOR AXIS 

= X-L '“CATION ALONG THE LINE OF LOADING CN THE WING 
= TOP SUCS'-VL MOMENT ARM OF THE 4-TH^ SECTION 
= SFA’JWI RE LENGTH CF EACH PANEL 
= C!j\Dr"NS30 CIRCULATION MATRIX 


-P-- 

O 


373 


C 

C 


pP r ( J > 


=Cf STANCE Of DYNAMIC PRFSSU4E OFF FPOM THE MAJOR AXIS 


379 

390 

3«1 

382 

3«3 

334 

335 
3«6 


3B7 

399 

3«9 

390 

391 

392 


393 

394 
393 

396 

397 

3<3ft 

349 

400 

401 

402 

403 

404 

405 

406 

407 

408 

409 

410 

411 

412 
4L3 

414 

415 

416 

417 
AT 3 

419 

420 

421 

422 


423 

424 

425 

426 
4?7 


IMPLICIT (A-H.0-71 

C DIM FUSION ETAlHPl) 

DIMENSION £TA(9) ,0PC(8) ♦FL0t48t 

L'!- w ?MSinw ARM! 31 tPCB.10) , XAX(8)iXL0l8)»Y|8) 

PIMrNSICN CLAC 8 i tOlfOEL CP<48,10» 

CpyupM /«„FA 2 / TfETAI 

COUPON M C [A 4 / R£ 0 INV{ 43 , 601 , A( 24 . 30 > 
r r M MQN ft Rf 45 / QELA< fa, lo J t FL C 6.101 

c 

2 Fr°MflT l«Fl:). 6 > 

8 FP°«AT I t 4X,* ACTUAL MOMENT ARM OF EACH COLUMN (ALONG THE SPANI ■ 
l I 

C 

THETA 2= 0.703 A? 

PPPTLOO.a 
9 r.pTT«= 10.0050 
T 1 9TF =10.0050 
TI "I* -5.900 
B*21.P-» 

e 

C CALTOLATr 01 STANCE OFF FROM THE MAJOR AXIS O p EACH STATION 
MRTsM +1 

="♦7 

nr 30 .Jsi.MPl 
CT4< jj-R=r( j-n t'A 
30 COM 7 IN' f 5 
0040 J=l,* 

Y ( J) - 0 . 50 * K TA CJ 1 +ET 4 C J + ll J 

Y ix 1. 1) =R EC Ttr *0. 5750 * < F COTTE-RPOTlF >*V l J»*D?4Nt THE TAll 
X in J> ap COT LC +0.250 * ( R COTTE-ROOTLE ) +YI J 1 *0TAN t The TA2I 
-*=!*MJ1 = X SXIJI-X LG<J» 

40 r7*.-INUE - * 

I >»N=M*N 

00 43 1=1.. MM 
OP 43 J= 1 .MD 2 
0: L CP( I*. >1=0.0 
OP 42 »=!,*] 
t=<v-l)avp? 

rr l CPIl.JleDFL CPU ..U-RFGINVII ,J*U 

42 CP*" i't'JF 

43 C r *TlNU7 

or 45 i = i.m 

r\r 45 J = 

Cl MI ,Jt=0.0 
?p 44 Kil. M 
L~ ( v — T > *V 

CL-tl. ) ) =CtM I tJl OcL CPII*l, J)/ 8 FTA*DELAIK, I + 2I 

44 CC'MIVjr 

r (T.J1=CLA( I. J )* SO/ 144.0 

45 CENT I NOE 
C 

: =:*‘p TW 7 EXACT LOCATION 09 DYNAMIC ORFSSURE 

‘■*:F5 = vrp*N 
?1r ’ 60 1=1. MN 

•LOf 1 1 = 0.0 

PC feO J=l t MN=S 

FL-m = 9L0m-RE0tNVtI ,J» 



■f 

H 


4? ft 
420 

430 

431 

432 

433 

434 

435 
435 
437 
43 B 

439 

440 

441 


44? 

443 


444 


445 

446 

44 7 
44fl 

449 

450 

451 
45? 

453 

454 

455 

456 

45 7 
4 SB 


60 CONTINUE 
DO 65 .1*1, M 
REP1=0.0 
. F.FP?=0.0 
On 62 I*l,N 
K=II-U*pi+J 

FL n l 1 1* FL0<5 1*nEL All t J+2 J/BSTA 
FLOCK )= FLDl K )* 50 2144.0 
p(spl=n?Pl + FLD(PI 
RrP2*R r ‘»»2-FLU »J*-2>*FLDiKt 
62 ‘ CONTINUE 

Ot>Cl JI*P.EP2/P3P1 
rPC(J)=QPC< J)*1.20 
65 CONTINUE 
C 

C PRINT OUT TM9 K0«5NT Af»M 

PR I NT R 

PRINT 2, HJPCC J) , J=1,H> 

C 

c COMPARE TO THE OUATFP CHORO NOMFNT AP.B 
PflNT ? , lfi6'*l J> ,J=1,B) 

c 

C F0<5N COEFFICIENT matrix F0» EXTERNAL LOAD 
00 50 1 

or so j = i,**stap 
iU.Jl-o.O 
50 CCNTTN'JE 

0^ *,9 t=l,N 
JfM JI=C D Cl I I 
T 3=1*3 

0<~ 55 J = 1.N£»2 
J3= J*3 

All ?-2 » J? I =r ( I , I) 

) = F( I, I) *000 5 1ALPHA) 

55 C : I *1>J - 

t Jt'.C.'j 

r'jn 

r m <»;rn iit-fuff ■ **«,= r * *■ * t* **■****=«:** * **» » * **»»**>= tr ***** *w ******* 




459 

460 
4M. 


SLPRDUTINF INVPTt A, N, SCALE) 

ILLICIT PFiL*a l 7 -H,T— £ ) 

ol'-'FMMnN AIN, M , TNOcX 190,2 ), IPVOTIROJ , P I VOT t 90 ) 
CALL *-paS"T( 207,256, 0,1) 



r 

5C1LF 3CWN 1 

462 


nr 1 1=1. N 

46? 


Dr l J=l,N 

464 


1 1IT.J)=A(1,J 

465 


067=1.0 

466 


7 1^=0. 0 

467 


67 0 r T = ' f : c 

465 


' 17 J = 1,N 

469 


17 ?CV r 7(J)=Q 

4 70 


PC 135 1=1, N 


C 

-“U°9T\ir i? 

4n 


*• =. * re -1 

4 7? 


« i=l,.V 

47 - 


"ir^o'ui- 

4 7 4 


" ? 3 5 * 1 ,'; 

475 


: s 1 ■ =>vr r <« 1 - 


.’3 


A 76 


43 

IFtfiAPS (T ) -DAB S f A f K ) 1 ) 83t23,?3 




*77 

.. 

03 

TPr W= j 




4 78 



iroi S K 




479 



T = UJ,K) 




480 


?3 

r. OMT! MUG 




401 


9 

c^mtinje 




48 2 



1 PVHT UCOL t-iovotf ICOL l + l 





C 


FOLLOW I MG 15 STATEMENTS PUT °IVOL ELF«FNT C!N DIAGONAL 




48? 



irMPOW-tCOL 1 73, 109 1 73 




434 


73 

r)FT=-PirT 




4P5 



nr 12 L =1 **i 




4«6 



T-4f TCriW,L) 




48 7 



AltBflW,! I=.M ICCL.U 




488 


12 

ACICri, U=T 




489 


109 

liF'OSXI ! ,l)=|P9U 




490 



It rr*< r , 2 ) =TCf}l. 


* 


441 



p Ivor (n=4ucai, icon 




492 



n''T=o?=T=oivuTC ! \ 



„ 


c 


F-\t •'•-‘IMG (> 5T/ rr^'.^TS T" 01 VI 05 POVOT ROW BY PIVOT ELEMENT 




49? 



if Kol, icol »*riM5 




494 



00 ?05 L = l,N 




495 


205 

A < ICPL . L > =' ( ICCLtL)/PI VOT( U 





c 


FOLLOWTNT 10 STA T EVENTS Ti? REDUCE NOM-PIvnT ROWS 




496 



ni P5 L!*1,U 




49 7 



IFfLI-ICCO 21,125.21 




499 


•21 

T = Ml T, ICUi 




499 



fiCL!.IC0L»=7tRC 




500 



r)r no l =1 .M 




501 


89 

A(ll, LI -Ml liLHAl ICOL ,L»*T 




50? 


135 

c nfis r i f ':'j ~ 





c 


^''Lir^rtT 11 STiT EVENTS T r IfJT'RCHAMGF COLUMNS 




503 


22 2 

C'C 3 1 = 1 




504 



L=N-t + 1 




505 



(U, I i-INO.£AU.2 5 ) 15,3.19 




506 


19 

j<n* = l\> x(i .,ii 




50 7 

l 


J'--'L=!'.r< , -Xf ,2) 




505 



S 4 - ifcj , \ 




509 







510 



A , J93W1 -i {K, JC0L1 




511 



•' IK, J'CLi -T 




512 


544 

C ■**;TI 




513 


3 

rr\: T Mj> 




514 

f* 

81 

C n *»T t M J E 





c 


SCALE 3*55 IMVRSE OF M5 TP 1 X-A-~v„__ i 




515 



Or 2 1 = 1 , M ‘ -* - 




516 



cr ? J=L,‘ 




517 


?. 

U i . J> = Atr,J?/SCAL£ 




519 



RETURN 





c 

81 

o?-u s N 




519 



r NO 





// 

1 -i 7 

A 





I\S'JPFIC IrCT MT”CPY t0 ASS ION* ARRAY STCFAGF* JOB A RAN CONED 




cnci UStjE 


TF J EC T COOS* 26312 OYT'StASSAY Af£A= 82752 RY TFS, TOTAL APE.* AVAILABLE* 

104544 

8YTFS 

DIAGNOSTICS 


O'- 1 , .ju^FQ OF WARNINGS* 1. NUMR.EP CF 

EXTFNSICNS= 

6 

CO IO 

Tt5F = 

6.69 S'-C •CjTtFiN Tl; 8 c = 0.00 SEC, T 6 M U7K*TFIV - VER l LEV 3 

JANUARY 

1572 

OAT£= 74/217 


-f . 

to 


44 


APPENDIX B 
T abulat ion, of Wind Tunnel Data 



= 423.06 sq. in., C = 

= 8.36 in. 


I) 50 




LEADING EDGE 

STATION NUMBER 

VERTICAL DEFLECTION 

a - 2,13 

ON WING 

5 le j 

fi TE 

; c T = 0.1683 

Jj 

3 

0.007 

-0.015 

C D = 0.103T 

4 

0.036 

-0.015 

C PM = -0.168U 

5 

I 

0.083 

, J 

0.042 

I 

TRAILING EDGE | 

6 

0.176 

■ 

. 

0.196 

' a = 2.13 : 

T 

0.301 

0.239 

; c * 0.1679 

) 

i 

8 

0.45 

0.301 

| c D = 0.1016 

9 

0.584 

0 . 486 

0 pM = -0.1635 | 

10 

0.731 

0.682 
l 





! LEADING EDGE 

STATION NUMBER 

VERTICAL DEFLECTION 

{ 

■ a ~ 5 . 37 

ON WING 

6 le 

6 te 

-c T = 0.1+61 

t. h 

3 

-0,029 

0.028 

= 0.1215 

» U 

4 

0.064 

0.079 

•C PM » -0.3301 

5 

0,21 

! 

1 0.163 

1 

1 

1 

: TRAILING EDGE 

6 

| 0.442 

1 

j 

‘J 0.379 

1 

i 

; a - 5.3T 

7 

0.752 

0.662 

i 

;C T = 0.1+6T5 

r i; 

8 

| 1.174 

\ 0.985 

,1 

: ;c = 0.121+2 

1 u 

9 

j 1.617 

1 1.1*56 

* 

j c pM - -0.331+1 


I 2.096 

! 

| 2 . 00 ^ 
1 


1*5 


LEAOTG EDGE 

STATION NUMBER 

VERTICAL DEFLECTION 

a = 8.6l 

ON WING 

5 le 

* 

TE 

C L = 0.758 

3 

0.021 

0.058 

C D = 0.11+95 

4 

| 0.19 

0.163 

C PM = -0.1*591 

5 

0.442 ; 

0.395 

TRAILING EDGE 

I 6 

0.849 

I 0.71 

a « 8. 6 

7 

l.4o6 

\ 1.182 

C T = 0.7521 

L ! 

8 

1 

2.089 

i 1.724 

c D - 0.1505 

9 

j 2.905 

j 2.491 

C PM - -0.1*552 

10 

j 3.742 

| 3.258 


II) a - 80 
00 


’LEADING EDGE 

*. 

station number 

VERTICAL DEFLECTION | 

: a = 2.13 

l 

ON WING 

S LE 

6 ^ S 
TE i 

'C_ = 0.1676 

' L 

3 

0.021 

-0.22 1 

1 

fc D “ 0.1059 

4 

o.o 64 

0.035 ■ 

jCpM = -0.1671 

! 

5 

0.132 

0.098 . 

i 

■ TRAILING EDGE 

J 

6 

0.282 

0.203 j 

! a = 2.13 

7 

0.463 

0.324 j 

c L = 0.167 

8 

0.718 

0.457 j 

•C^ = 0.1033 
D 

9 

0.908 

0.747 i 

^PM ” 

L - • - 

1.175 . 

1.019 ’’ 













